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Abstract 

We discuss two methods by which high porosity silica aerogels can be engineered to exhibit global 
anisotropy. First, anisotropy can be introduced with axial strain. In addition, intrinsic anisotropy 
can result during growth and drying stages and, suitably controlled, it can be correlated with 
preferential radial shrinkage in cylindrical samples. We have performed small angle X-ray scattering 
(SAXS) to characterize these two types of anisotropy. We show that global anisotropy originating 
from either strain or shrinkage leads to optical birefringence and that optical cross-polarization 
studies are a useful characterization of the uniformity of the imposed global anisotropy. 

PACS numbers: 61.05.cf, 67.30.hm, 78.20.Fm, 82.70.Gg 



I. INTRODUCTION 



High porosity silica aerogel is a transparent, low density, material consisting of aggre- 
gates of silica particles, often described as strands, with a diameter of ~ 30 — 50 A and 
average separation, or correlation length, £ = 100 — 1000 ApQ. The microstructure of the 
aerogel is inherently inhomogeneous on the scale of the correlation length, which we take 
to be the upper cut-off of a fractal distribution of the silica particles. Fractal distributions 
follow naturally from the physical process of diffusion limited aggregation that occurs during 
gel formation and have been characterized with small angle scattering experiments (visible 
light [2J, X-ray[3], and neutrons jl]) as well as having been simulated numerically|5j|6j. In the 
present work we are particularly interested in long length scale inhomogeneity, specifically 
the contolled introduction of global anisotropy into silica aerogels and its characterization. 
In the following we will use the term isotropic to refer to an aerogel having no preferred 
direction on average over macroscopic specimens that are typically of mm size or larger. 

Silica aerogels with porosities up to ~ 99% are of interest in a wide variety of scientific 
and engineering disciplines including high-energy physics, astrophysics, material science, 
condensed matter physics and chemical engineering. Aerogels find application in particle 
physics to detect Cerenkov radiation [7j and have been used by NASA to collect particulate 
matter from comet 81P/Wild 2 [8] [9]. Silica aerogels have seen widespread application in 
materials science and because of their fractally correlated structure [3] have basic scientific 
significance. Aerogels are used in the manufacture of ceramic materials [TU]. They are impor- 
tant in the study of critical phenomena and exploration of impurity effects in quantum fluids 
such as 4 He pi] and 3 Hep2] as well as to understand the influence of impurities on liquid 
crystal phase transitions [13] [1.4]. In addition, hybrid materials based on silica aerogels can 
be engineered with metal, or metal oxide, particles incorporated into the aerogel framework 
for the purpose of producing heterogeneous catalysts [13]. 

Global anisotropy of the aerogel is an important consideration in many of these systems. 
In general, structural characterization is necessary to understand the behavior of fluids in 
porous media such as aerogels. This is also true for anisotropic structures which can lead to 
anisotropic diffusion and flow and therefore be of consequence in many of the aforementioned 
applications. For example, the ordered phases of liquid crystals in aerogel are inherently 
anisotropic [TB"] and we can expect that they will be modified by global anisotropy [TB]. The 
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physical and optical properties of aerogels are modified by anisotropy, as we discuss in 
this work, and this may be important in some applications such as transparent thermal 
insulation [T7j. More generally, anisotropic aerogels are members of a larger group of fasci- 
nating anisotropic porous media including sedimentary rock[T8], bone|19j. and cementitious 
materials [20j. 

Our motivation for the present study is linked to our interest in a basic physics problem 
concerning impurity effects on superfluid 3 He. It is thought that anisotropy will play a role in 
determining which phases of superfluid 3 He in aerogel are energetically stable |2H [22] and that 
anisotropy will influence the formation of textures of the superfluid order parameter. Here 
we report that high porosity silica aerogels can be engineered to exhibit global anisotropy 
and we discuss small angle X-ray scattering (SAXS) and optical polarization methods for 
characterization of this anisotropy. 

II. SAMPLES 

The aerogel samples discussed in this work were grown at Northwestern University via the 
"one-step" sol-gel methodj23]. A silicon oxide precursor, tetramethyl orthosilicate (TMOS), 
is dissolved in methanol and hydrolyzed. Ammonia is used as a catalyst. 

Si(CH 2 OH) 4 + 2H 2 -> Si0 2 + 4CH 3 OH 

This reaction results in a wet gel, known as alcogel. The alcogel is then supercritically dried 
following a method that is based on the rapid supercritical extraction process (RSCE) of 
J.F. Poco et al. [21]. 

Specifically, the sol is poured into a cylindrical stainless steel chamber having an inner 
diameter of 5.08 cm and a length of 4.92 cm and containing glass tubes of various sizes as 
molds. The chamber is completely filled with the sol and sealed using two stainless steel 
lids each having a thickness of 1.27 cm. Each lid is bolted in place using six 8-32 screws. 
The alcogel that forms is typically aged at room temperature for approximately 3 days. The 
chamber containing the alcogel is then loaded into an autoclave having an open volume of 
0.46 L and containing ~ 50 mL of methanol. The autoclave is loaded into an oven and 
heated through several manual adjustments and, as the temperature rises, the pressure in 
the alcogel chamber forces the lids to open slightly, allowing fluid to escape into the open 
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volume of the autoclave. Once the critical point of methanol (T c = 239. 5°C, p c = 80.81 
bar) is exceeded the autoclave is emptied by slowly opening a high pressure valve. This 
drying typically begins from a starting temperature and pressure of ~ 270°C and ~ 100 bar. 
The system is depressurized ~ 8 — 12 hours and then allowed to cool to room temperature 
overnight. Fig. [T] depicts the typical oven temperature and autoclave pressure profiles during 
aging and drying. With this method we have successfully produced silica aerogel ranging in 
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FIG. 1: (Color online) Oven temperature and autoclave pressure profile during aging and drying. 



porosity from ~ 94% — 99%. A sample's porosity is determined from its weight and volume 
after supercritical drying to an accuracy of 0.1% for a 98% aerogel. 

For the SAXS and optical birefringence measurements reported here, cylindrical samples 
with porosities ranging from ~ 97% — 98% were grown in glass tubes with ~ 8 mm inside 
diameter. By controlling the catalyst concentration and the supercritical drying procedure 
we could grow samples with radial shrinkage in the range 0% to 20%. From these samples 
we chose those with no significant shrinkage on-axis. The samples were cut to ~ 0.75 cm 
with aspect ratios of ~ 1. 



III. SAXS EXPERIMENTS 



The scattering experiments were carried out at Sector 8 of the Advanced Photon Source 
(APS) at Argonne National Laboratory, using a photon energy of 7.5 keV, corresponding to 
a wavelength, A = 1.69 A. The X-ray beam, with a spot size of ~ 100 microns, was incident 
on the the center of a sample, which was oriented with its cylinder axis perpendicular to 
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the beam. The experiments discussed below can be grouped into two categories. First, 
to study the effect of axial strain on the aerogel, scattering experiments were carried out 
on a nominally isotropic sample as well as on a sample exhibiting 13.8% radial shrinkage. 
They were subjected to increasing strain in the range 0% to ~ 32%. Secondly, in order to 
systematically explore the relationship between radial shrinkage and X-ray scattering, we 
conducted SAXS on a series of samples exhibiting radial shrinkage in the range 0% to 20%. 

In each experiment we recorded the scattered X-ray intensity, I(q), as a function of the 
momentum transfer, q, using a CCD camera. The momentum transfer is defined as q = 
(47r/A) sin(#/2), where 9 is the angle between the incident and scattered X-ray wavevectors. 
For a given image, we binned I(q) in ~10° increments of the azimuthal angle, 0, defined in 
the plane of the CCD camera (refer to Fig. [2]). The cylinder axis was defined to be = 90°. 
Fig. [2] is a cartoon depicting the geometry of a typical SAXS measurement. 



CCD Camera 




FIG. 2: (Color online) Sketch of the experimental geometry. Note that = 90° is parallel to the 
cylinder axis. 



IV. SAXS RESULTS AND DISCUSSION I: AXIAL STRAIN AND ANISOTROPY 

The dependence of I(q, 0) on provides information about the magnitude and direction 
of anisotropy in the aerogel. This is depicted in Fig. [3j which presents images of the scattered 
X-ray intensity from two aerogel samples. The sample labeled (a) is of a nominally isotropic 
and unstrained aerogel, while (a') depicts the same sample, but strained by 12.7% along the 
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cylinder axis. The image labeled (b) is of a sample exhibiting 13.8% radial shrinkage. In 
what follows, the labeling scheme presented in Fig.[3]will be used when referring to these two 
samples. The porosity of sample (a) was 98% and that of (b) 97.6%. Fig. [3] qualitatively 

Min Max 



(a) (a') (b) 

FIG. 3: (Color online) False color scattered X-ray intensities for (a) a 98% nominally isotropic 
aerogel, (a') the same sample unaxially strained by 12.7%, and (b) a 97.6% aerogel exhibiting 
13.8% radial shrinkage. In each case the sample is oriented with its cylinder axis vertical relative 
to the images. The dark rectangle is a beamstop inserted to protect the CCD camera from the 
unscattered transmitted beam. 

demonstrates the existence of two distinct types of anisotropy. First, comparing (a) and 
(a') shows that axial strain can be used to introduce anisotropy into a nominally isotropic 
aerogel. In addition, (b) indicates the presence of intrinsic anisotropy in samples exhibiting 
radial shrinkage. Moreover, this intrinsic anisotropy is ~ 90° out of phase with that induced 
by axial strain. 

We analyzed the scattering curves by fitting I(q, 0) with the following phenomeno logical 
scattering function, 

CP d C1 4- a 2 ? 2 ) 1 / 2 
%0 = (1 + ^/4 J sinP-Dtan-^OL (1) 

where C, d, and £ are fit parameters which depend on 0. Eq. 1 is a modified version of 
a structure factor for fractally correlated materials having an upper length scale cut-off. 
The (/-dependence of Eq. 1 is stronger than the structure factor derived by Freltoft et. 
al. [25]. This sharper g-dependence was chosen to better match our experimental data. The 
parameter d is the fractal dimension of the aerogel, and C is an overall scale factor. The 
length scale £ is associated with the aerogel correlation length, i.e. the upper length scales 
at which the aerogel ceases to be fractally correlated and, for future reference, we note that 
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the axial (radial) correlation length corresponds to = 90° (0 = 0°). Note that in the limit 
g£ 3> 1, Eq. 1 is proportional to q~ d as expected [23] for the fractal regime. We performed fits 
to our data using Eq. 1 and obtained an excellent representation as can be seen in a typical 
case in Fig. [4] for sample (a') (0 = 5.3°, = 88.5°). Fig. [4] shows that when the aerogel is 
compressed along the cylinder axis the axial correlation length is compressed relative to the 
correlation length in the radial direction. 
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FIG. 4: (Color online) I(q,<j>) fit with Eq. 1 (0 = 5.3°, = 88.5°) for the sample (a'). The curves 
have been offset vertically for clarity, otherwise the data points would coincide at high q. 

To quantify the magnitude and direction of the anisotropy in a given sample, we plot the 
correlation length £ versus 0. For all of the samples we found the following 0-dependence, 

£(0)=£o + 6cos(20), (2) 

where £o and £i are fit parameters. Note, the average fractal dimension ranged from 1.55-1.80 
over the samples discussed in this work. For any given sample, the fractal dimension, d, was 
found to be a much weaker function of 0, as compared to £(0). The dependences of d and C 
on are not important to our analysis and do not pertain to our discussion. The parameter 
£o represents the average correlation length of the aerogel and was found to be of the order 
of ~200 A. The amplitude, £i, is a measure of the magnitude of the anisotropy. In Fig. [H] 
we present this 0-dependence for the samples (a), (a'), and (b). The sinusoidal dependence 
of £(0) demonstrates that anisotropy from axial strain, or radial shrinkage, exists on the 
length scale of the aerogel correlation length, £. Moreover, comparing the traces for samples 
(a') and (b) demonstrates that the axial strain induced anisotropy has the opposite sign of 
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FIG. 5: (Color online) £-£o as a function of (f> for the samples (a), (a'), and (b). The solid curves 
are fits using Eq. 2. The amplitude of the sinusoid is a measure of the magnitude of anisotropy. 

that from anisotropy produced by radial shrinkage. Fig. [5] also indicates that the nominally 
isotropic sample, i.e. sample (a), also possesses a small amount of intrinsic anisotropy, which 
has a magnitude d ~ 9 A. 

To further investigate the relationship between axial strain and global anisotropy we plot 
£i versus applied axial strain. The results for samples (a) and (b) are presented in Fig. [6j 
The increase of £i with increasing strain demonstrates that anisotropy can be introduced 




60 


i 


i i i i i i 


40 
20 





. . ■ ■ * 

* 


♦ 


■ 


-20 




■ 

♦ a (nominally isotropic) 


-40 


■ 


■ b (radial shrinkage = 13.8%) " 


-60 


i 


i i i i i i 



5 10 15 20 25 30 35 
Strain (%) 

FIG. 6: (Color online) £i as a function of axial strain applied to samples (a) and (b). The vertical 
offset between the two data sets reflects the difference in magnitude of the intrinsic anisotropies of 
the two samples. 
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into the aerogel systematically. Moreover, axial strain can be used to tune between the 
two types of anisotropy and implies that strain can be used to compensate for intrinsic 
anisotropy. The negative values of £1 are an indication of the overall phase in the second 
term in Eq. 2. 

We also plot £o versus strain to see how the average correlation length of a sample is 
affected. The results for samples (a) and (b) are presented in Fig. [7| We find that £o 
decreases monotonically with applied strain for both samples. Previous SAXS studies on 
isostatically strained aerogels also indicate a decrease in the correlation length with applied 
strain! 
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FIG. 7: (Color online) £o as a function of axial strain applied to samples (a) and (b). 



Fig. [6] and Fig. [7] seem to suggest that the primary effect of strain is to reduce the axial 
correlation length while only weakly increasing the radial correlation length. The correlation 
length parallel (perpendicular) to the cylinder axis can be obtained from Eq. 2 by setting 
= 90° (0 = 0°) and we define £y = £o — £i = £o + £i) as the axial (radial) correlation 
length. In Fig. M we plot £y and £j_ divided by their values at zero strain, £|| and £j_o, 
versus strain applied to samples (a) (uppper panel) and (b) (lower panel). Fig. [8] clearly 
demonstrates that axial strain primarily compresses i.e. the correlation length parallel 
to the strain axis, while only modestly increasing £j_ . 

Fig. [8] also allows us to determine the extent to which the macroscopic strains are trans- 
mitted to microscopic length scales, i.e. length scales of the order of £. For small values of 
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FIG. 8: (Color online) Normalized axial (radial) correlation length, £||/£||o (£l/£lo)> versus strain 
for sample (a) (upper panel) and (b) (lower panel). The solid lines are linear fits to the low strain 
portion of C|l/C|lo an d £±/£_i_o- The dashed line corresponds to a = 1 in Eq. 3 and 4. 



strain the data in Fig. [8] indicate, 



1 + aeu 



fie* 



(3) 
(4) 



where ey = A///o is the axial strain on a cylindrical aerogel intially having a length of Iq and 
compressed by A/. The constants a and /3 describe the transmission of macroscopic strain 
down to the level of the correlation length and are determined from linear fits of £||/£||o an d 
£_l/£_lo near zero strain. These fits are presented in Fig. [8] as solid lines for samples (a) 
(upper panel) and (b) (lower panel). In Eq. 3 and 4 we use the convention that e\\ < for 
compression. Homogeneous elastic theory (HET) suggests that a = 1 and (5 = v, where v 
is the macroscopic Poisson ratio. Note, a = 1 implies that the macroscopic strain scales 
directly down to the level of the aerogel correlation length and is represented by the dashed 
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lines in Fig. [8j If we had found a < 1 it would suggest that strain is not distributed to 
microscopic length scales, while a > 1 would imply an enhancement at microscopic length 
scales. The dependence of £||/£||o on strain in Fig. [8] unambiguously indicates a value of a > 1, 
implying that the aerogel cannot be accurately described by HET. We find a = 2.5 ±0.1 
for sample (a) and a = 2.9 ±0.1 for sample (b). However, from our measurements we find 
that the ratio /3/a « v. Specifically f3/a = 0.20 ±0.01 for sample (a) and (3/a = 0.18 ±0.02 
for sample (b). These values are qualitatively consistent with typical values of v for silica 
aerogels[27j. In addition, we have measured the macroscopic Poisson ratio for a sample 
similar to (a) using optical techniques (see section VII) and find v = 0.30 ± 0.05. 

We conclude that aerogel is not a homogenous elastic medium and that axial strain 
induces larger effects on the correlation length scale than would be expected for such a 
medium. Finally, a Poisson-like behavior describing the microscopic displacements appears 
to hold. 



V. SAXS RESULTS AND DISCUSSION II: INTRINSIC ANISOTROPY AND RA- 
DIAL SHRINKAGE 

To explore the possibility of systematically introducing anisotropy different from that 
produced by strain, we have also performed SAXS on a series of aerogels which exhibited 
preferential radial shrinkage after supercritical drying. Following the same analysis as de- 
tailed in the previous section, we studied the evolution of the the anisotropy amplitude, 
£i, as a function of radial shrinkage. The results are presented in Fig. [9j Specifically, the 
absolute value of £i appears to increase systematically with increasing shrinkage, indicating 
a correlation between radial shrinkage and anisotropy. The larger data scatter in Fig. 7 as 
compared with Fig. 5, reflects that our measurement and control of radial shrinkage is not 
as uniform and systematic as for imposed strain on a single sample. However, the small 
anisotropy of the samples having small shrinkage is prominent. It is noteworthy that the 
magnitude of the effect on |^| is comparable with that observed for strained samples. 

In addition, we plot the average correlation length, £ , as a function of radial shrinkage in 



Fig. 10 For the samples that exhibited radial shrinkage, £o was found to decrease by as much 
as ~40%. This effect is of similar magnitude as that observed for strained samples with the 
exception of the samples that did not show any shrinkage and which are not consistent with 
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FIG. 9: (Color online) Anisotropy amplitude, £i, as a function of radial shrinkage. The dashed line 
is a guide to the eye. Note the relatively small magnitude of intrinsic anisotropy in the samples 
that did not exhibit shrinkage. 
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FIG. 10: (Color online) Average correlation length, as a function of radial shrinkage. The 
dashed line is a guide to the eye. Note, the behavior of the samples exhibiting no shrinkage is 
inconsistent with the general trend presented here, for which we have no explanation. 



the the general trend in Fig. 10 , for reasons we do not understand. 



VI. OPTICAL BIREFRINGENCE EXPERIMENTS 



When a transparent material possesses an anisotropic dielectric constant it will exhibit 
birefringence. Illuminating an optically birefringent material with white light will result in 
transmission of two components with orthogonal polarizations, parallel and perpendicular 
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to the "optical axis". By placing such a material between crossed polarizers it is possible 
to determine if there is a well defined optical axis and how it is distributed throughout the 
sample [25] . Measurements of optical birefringence are well-known for characterization of 
liquid crystals [29] and many other materials [30] . 

We have performed optical birefringence experiments on high porosity aerogels. Cylindri- 
cal aerogel samples were placed between two crossed polarizers and illuminated with diffuse 
white light. The cylinder axis was oriented vertically, a geometry similar to that shown in 
Fig. [I] A diffuser was placed between the light source and the first polarizer and images 
were recorded with a digital camera located after the second polarizer. 

To study the effect of axial strain on a nominally isotropic aerogel, optical birefringence 
experiments were conducted as the sample was subjected to increasing strain in the range 
0% to 18.6%. In addition, we have investigated samples with radial shrinkage and have 
confirmed that optical birefringence can be used to detect global anisotropy in this situation 
as well. In both cases, the optical axis is primarily oriented along the cylinder axis. 



VII. OPTICAL BIREFRINGENCE RESULTS AND DISCUSSION 



Fig. 11 presents optical birefringence of a 98% porosity aerogel as it was axially strained 
to 18.6% in 2.3% steps. The crossed polarizers were oriented with one polarizer at 45° and 
the other at 135° relative to the vertical cylinder axis. This sample did not exhibit any radial 



shrinkage after supercritical drying. Panel 1 of Fig. 11 indicates that the unstrained gel is 
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FIG. 11: (Color online) Optical birefringence of a 98% porosity aerogel, nominally isotropic before 
it was subjected to increasing axial strain. The strain increases from left to right: 1 (unstrained), 
2 (2.3% strain), 3 (4.7%), 4 (7.0%), 5 (9.3%), 6 (11.6%), 7 (14.0%), 8 (16.3%), 9 (18.6%). 



homogeneously isotropic since it does not impose any preferred direction of polarization 
when placed between crossed polarizers for all orientations of the polarizers with respect 
to the cylinder axis. We have also found no transmission of light propagating down the 



13 



cylinder axis for this isotropic aerogel placed between crossed polarizers. Panels 2-9 of 
Fig. [TTJ demonstrate that strain can be used to effectively convert this nominally isotropic 
aerogel into a polarizer. We attribute this effect to optical birefringence. The polarizing 
effect increases with increasing strain and is the result of structural anisotropies on optical 
length scales. Furthermore, it is evident from the approximately uniform intensity of the 
image that this anisotropy is a global property of the entire sample and is homogeneously 



distributed. From the images presented in Fig. [TT] we can measure the relative intensity of 
light passing through the sample averaged over a central region of the image of the aerogel as 



a function of strain. These results are presented in Fig. 12 and indicate that the transmitted 



intensity increases linearly with increasing strain up to the largest value of strain, i.e. 18.6%. 
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FIG. 12: (Color online) Relative intensity averaged over a central region of the image of the sample 



versus strain for the sample presented in Fig. 11 



By rotating both polarizers, keeping them crossed, we can determine the direction of the 



anisotropy axis. Fig. 13 presents such a rotation sequence for the sample at it's maximum 
compression, i.e. 18.6%. The fact that the intensity maxima (minima) are seen when the 
polarizers are oriented at 45° and 135° (90° and 180°) is consistent with the anisotropy axis 
being oriented along the cylinder (strain), axis. 

To complement our SAXS measurements, we also used optical birefringence to investigate 
possible intrinsic anisotropy in a sample exhibiting radial shrinkage. As before, the polarizers 
were intially oriented with one at 45° and the other at 135° relative to the cylinder axis. They 
were then rotated together, keeping them crossed. We have found that aerogels exhibiting 
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FIG. 13: (Color online) Optical birefringence of a 98% porosity aerogel strained by 18.6%. The 
labels are associated with the rotation of the polarizers relative to the cylinder axis: 1 (45°, 135°), 
2 (60°, 150°), 3 (75°, 165°), 4 (90°, 180°), 5 (105°, 195°), 6 (120°, 210°), 7 (135°, 225°). 

radial shrinkage polarize transmitted light and this can be observed with cross polarizers, 
indicating the existence of a well defined optical axis. Again we attribute this effect to 



optical birefringence. In Fig. 14 we present the results for an aerogel cylinder exhibiting 
12.7% radial shrinkage showing that axial intrinsic anisotropy is present in this sample. It 
is noteworthy that this intrinsic anisotropy seems to be less uniform throughout the sample 



than that induced by strain as displayed in panels 2, 3, 5, and 6 of Fig. 14 
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FIG. 14: (Color online) Optical birefringence of a 98% porosity aerogel exhibiting 12.7% radial 
shrinkage. The labels are associated with the rotation of the polarizers relative to the cylinder axis: 
1 (45°, 135°), 2 (60°, 150°), 3 (75°, 165°), 4 (90°, 180°),5 (105°, 195°), 6 (120°, 210°), 7 (135°, 
225°). 



For the intrinsically anisotropic samples, i.e. those exhibiting radial shrinkage, we ob- 
served only a small amount of transmission through crossed polarizers for light propagating 



down the cylinder axis. This, along with the rotation series presented in Fig. 14, indicates 
that the cylinder axis is predominantly the optical axis in the case of samples exhibiting 
anisotropy due to radial shrinkage. 

Our SAXS measurements indicate that microscopic anisotropy is associated with strain 
or radial shrinkage developed during the sample growth and drying stages. Furthermore, 
with a SAXS spot size of ~ 100 microns it is possible to image the distribution of anisotropy 
by scanning the sample profile. We have found that tabletop optical polarization studies 
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are complementary and less time consuming. They give a clear picture of the principal 
directions of the anisotropy on optical length scales and have the advantage of displaying an 
image that shows the homogeneity of this anisotropy as it is distributed over the sample. 

VIII. CONCLUSIONS 

In conclusion, our measurements demonstrate that high porosity silica aerogels can be 
engineered with global structural anisotropy produced either by applied strain or during 
sample growth and drying. SAXS measurements indicate that this anisotropy is on the 
length scale of the aerogel correlation length, £, and optical birefringence measurements 
confirm that it is a global feature of the aerogel. These aerogels represent a novel type of 
anisotropic porous medium with possible implications for applications in a wide variety of 
physical systems. 
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